Remi Daigle

require(tidyverse)
require(sf)
require(leaflet)

I am working on 3 CHONe projects, because 1 is not enough.

Here’s the code I used to map out my study sites. For OHI, I have a shapefile, but for the others, I’m only creating a box around the coordinates.

# set standard projection
projll <- "+proj=longlat +datum=WGS84"
# set OHI project
OHI <- read_sf("shapefiles/regions_final.shp") %>% 
    st_simplify(preserveTopology=T,dTolerance=1000) %>%
    st_transform(projll) %>% 
    st_combine() %>% 
    st_sf(name="Remi Daigle",
         description="OHI regional assessment",
         CHONeProject="2.1.1",
         geometry=.) 
# set Atlantic LTRANS connectivity project
bcbbox <- c(-70.88808,39.89707,-47.28724,53.56291)
BC <- cbind(
            c(bcbbox[1],bcbbox[3],bcbbox[3],bcbbox[1],bcbbox[1]),
            c(bcbbox[2],bcbbox[2],bcbbox[4],bcbbox[4],bcbbox[2])
        ) %>% 
    list() %>%                  # put that matrix in a list
    st_polygon() %>%            # put that list in a simple features polygon
    st_sfc(crs = projll) %>%  # put that polygon into a simple features column and give projection info
    st_sf(name="Remi Daigle",
         description="Atlantic LTRANS connectivity",
         CHONeProject="1.1.1",
         geometry=.)%>%    # finally make it an sf object
    st_cast("MULTIPOLYGON")
# set Atlantic LTRANS connectivity project
Atlbbox <- c(-135.8003,45.3837,-122.0694,56.1647)
Atlantic <- cbind(
            c(Atlbbox[1],Atlbbox[3],Atlbbox[3],Atlbbox[1],Atlbbox[1]),
            c(Atlbbox[2],Atlbbox[2],Atlbbox[4],Atlbbox[4],Atlbbox[2])
        ) %>% 
    list() %>%                  # put that matrix in a list
    st_polygon() %>%            # put that list in a simple features polygon
    st_sfc(crs = projll) %>%  # put that polygon into a simple features column and give projection info
    st_sf(name="Remi Daigle",
         description="BC LTRANS connectivity",
         CHONeProject="1.1.1",
         geometry=.) %>%    # finally make it an sf object
    st_cast("MULTIPOLYGON")
# bind all my sites together
sites <- rbind(OHI,BC,Atlantic)

Finally, plot it!

cols <- RColorBrewer::brewer.pal(nrow(sites),"Dark2")
leaflet(sites) %>%
    addTiles() %>% 
    addPolygons(group = ~description,
                color = cols) %>%
     addLegend("topright", 
               colors = cols,
               labels = unique(sites$description))

NA